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We calculate the effect of gravitational wave (gw) back-reaction on realistic neutron stars (NS's) 
undergoing torque-free precession. By 'realistic' we mean that the NS is treated as a mostly-fluid 
body with an elastic crust, as opposed to a rigid body. We find that gw's damp NS wobble on 
a timescale re ~ 2 x 10^ yr [10~^/(A/d//o)]^(kHz/z/s)'*, where Us is the spin frequency and Aid 
is the piece of the NS's inertia tensor that "follows" the crust's principal axis (as opposed to its 
spin axis). We give two different derivations of this result: one based solely on energy and angular 
momentum balance, and another obtained by adding the Burke-Thorne radiation reaction force to 
the Newtonian equations of motion. This problem was treated long ago by Bertotti and Anile (1973), 
but their claimed result is wrong. When we convert from their notation to ours, we find that their 
Te is too short by a factor ~ 10^ for typical cases of interest, and even has the wrong sign for Aid 
negative. We show where their calculation went astray. 
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I. INTRODUCTION 

This paper calculates the effect of gravitational wave 
(gw) back-reaction on the torque-free precession, or wob- 
ble, of realistic, spinning neutron stars (NS's). By 'reahs- 
tic' we mean the NS is treated as a mostly-fluid body with 
an elastic crust, as opposed to a rigid body. (However 
we do not include any superfluid effects in our analysis.) 
Freely precessing neutron stars are a possible source for 
the laser interferometer gw detectors (LIGO, VIRGO and 
GEO under construction, TAMA already operational); it 
is the prospect of gravitational wave astronomy that mo- 
tivated our study. Also, the first clear observation of free 
precession in a pulsar signal was reported very recently 

, making this investigation all the more timely. 

The effect of gw back-reaction on wobbling, ax- 
isymmetric rigid bodies was first derived 27 years 
ago, in an impressively early calculation by Bertotti 
and Anile [||. They found (correctly) that for 
rigid bodies, gw backreaction damps wobble on a 
timescale (for small wobble angle 9) r^'®"^ = 1.8 x 
106yr[10-'^/(A///i)]2(kHz/i/s)4(1045gcm7/i), where 
is the spin frequency and AI = (/s — Ii) (with 

In the same paper, Bertotti and Anile went on 
to calculate the effect of gw back-reaction on wob- 
ble for the more realistic case of an elastic NS. When 
cast into our notation, their claimed gw timescale is 
5IiC^ /[2G{2TTVs)'^AlQAId] where AIq is the asymme- 
try in the moment of inertia due to centrifugal forces 
and Aid is the asymmetry due to some other mech- 
anism, such as strain in the solid crust. Taking 
A/o to be (roughly) the asymmetry expected for a 
rotating fluid according to A/o// « 0.3(z^s/kIIz)^, 
we would then have a damping time of merely 
0.6yr(kHz/z/s)^[10"7(A/d//)](10'*5gcmV/). Despite 
the fundamental beauty of this probem and its poten- 
tial astrophysical significance, their remarkable claim- 
that in realistic NS's, gw's damp wobble with amazing 
efficiency-was apparently little known. (A citation index 
search showed that Bertotti and Anile |^ had been ref- 
erenced by other authors only four times in the last 27 
years.) 

We will show that the Bertotti and Anile result for 
elastic NS's is very wrong, however. For typical cases of 
interest, their gw timescale rg is too short by a factor 
~ 10^ ! Moreover, their calculation even gives the wrong 
sign (exponential growth instead of damping) when Aid 
is negative. R In contrast, we find that gw's always act to 



damp the wobble in realistic NS's, just as for rigid bodies. 
While in Nature the typical case will be Aid positive. 
Aid < can also occur in principle. We call attention 
to this case not because it is common, but because it 
highlights how much our result differs from Bertotti and 
Anile and because, in fact, their implicit prediction 
of exponential wobble growth for this case provided our 
initial impetus to look more closely at this problem. 

The organization of this paper is as follows. In §2 we 
derive the gw damping timescale for rigid-body wobble, 
using the mass quadrupole expressions for the energy and 
angular momentum radiated to infinity. (This derivation 
is actually Exercise 16.13 in the textbook by Shapiro and 
Teukolsky |^). We give another derivation of Tg in §3, 
this time by adding the Burke-Thorne radiation reaction 
force directly to the Newtonian equations of motion. This 
latter approach was how Bertotti and Anile first cal- 
culated (correctly) the gw damping time for wobbling, 
rigid bodies. 

In §4 we review standard material on the torque-free 
precession of elastic bodies, in the absence of viscous 
terms or gw back-reaction. In §5 we derive the gw damp- 
ing timescale rg in the elastic case, using energy and an- 
gular momentum balance. In §6 we give a second deriva- 
tion of Tg in the elastic case, using the Burke-Thorne 
radiation reaction force to evolve the elastic body's free 
precession. This was also the strategy of Bertotti and 
Anile (I, and we show where they went wrong. Briefly, 
they did not realize that in addition to torquing the NS, 
the radiation reaction force also perturbs the NS's shape 
(in particular, its inertia tensor). When solving for the 
evolution of the wobble angle, we show that the "per- 
turbed shape" term in the equations of motion almost 
entirely cancels the gw torque term that they do include. 
(Of course, by definition there is no "perturbed shape" 
term in the rigid-body case, which is probably why they 
forgot this term when adapting that calculation to the 
elastic case.) In §7 we describe how to include the ef- 
fects of a fluid core in the radiation reaction calculation. 
Finally, in §8 we conclude by commenting briefly on the 
astrophysical implications of our result. 

We will work in cgs units. 

II. RADIATION REACTION FOR A RIGID 
BODY: ENERGY AND ANGULAR MOMENTUM 
BALANCE 

The derivation of the wobble damping rate for realistic 
NS's, using energy and angular momentum balance, is 



^Actually, Bertotti and Anile [g] never claim in words that 
they find unstable growth of the wobble angle when Aid < 0, 
but that is what is found if one just takes their formulae 
and converts from their notation to ours, as above. More- 
over we have repeated their (flawed) calculation, includ- 



ing their one crucial error, and seen that it does lead to 
a prediction of exponential wobble growth for Aid nega- 
tive. The conversion from their notation to ours is simply 
(Sil - (52/)(cos2 7 - i sin^ 7) ^ Ah and S2I -» A/n). 
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rather similar to the corresponding derivation for rigid 
bodies. Here we briefly review the solution to the rigid- 
body problem, as a warm-up for tackling the realistic 
case. 

Consider an axisymmetric rigid body with principal 
axes xi,X2, £3 and principal moments of inertia Ii = I2 =/= 
I3. Let the body have angular momentum J, misaligned 
from X3. Define the wobble angle 6* by J • X3 = JcosO. 
It is a standard result from classical mechanics that 
(in the absence of external torques) the body axis £3 
precesses around J with (inertial frame) precession fre- 
quency (j) = J/Ii, with 9 constant [Q. Together, the pair 
{0, (j)) completely specify the free precession (modulo a 
trivial constant of integration specifying at t = 0). We 
wish to calculate the evolution of these two parameters 
using the time-averaged fluxes [E, J). 

Straightforward application of the mass quadrupole 
formalism Q gives 



E = 



2G 

5^' 



'i)^{Mf sin^ 6l(cos^ 9 + IGsin^ 6*), (2.1) 
where A/ = I3 — Ii, and 



j^E/(f>. (2.2) 
It follows from differentiation of = J/ Ji that 

,;i = - — — — 05sin2 6»(16sin2 6'-hcos2 6l). (2.3) 

To calculate the rate of change of the wobble angle, re- 
arrange 
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to give 
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where Eq. (2.2) has been used. The energy of the body 
is simply its kinetic energy: 



and so 
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This gives 



5c5 h 



cos 9 sin e{lQsm^9 + cos^9). (2.9) 



We can construct timescales on which the spin-down and 
alignment occur: 
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he' 1 h 



2G 04 A/2 sin^ 6'(16sin^6'-f cos^l 



(2.10) 
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dt 



sin 9 2G ^4 A/2 cos2 61(16 sin^ 9 + cos^ 



(2.11) 



Radiation reaction causes both and sin^ to decrease, 
regardless of whether the body is oblate or prolate. Note 
that in the limit of small wobble angle the inertial pre- 
cession frequency remains almost constant ('''^'^"^ — > 00), 
while 9 decreases exponentially on the timescale 



rigid _ 



5c5 1 /i 



(2.12) 



Parameterising: 



rfsid = 1.8 X lOVr 



/ 10- 



VA///1 



lO''^ gcm2 



/i 



(2.13) 



In the limit of vanishingly smal l wo bble angle the par- 
tial derivative on the Ihs of Eq. ( |2.7| ) becomes what we 
conventionally call the 'spin frequency' of the body . 
Eq. ( ^.5[) then shows that 9 is proportional to the differ- 
ence between the inertial precession frequency 4> and the 
spin frequency fi. This difference remains finite as — > 
according to - O = {AI/Ii)n[l + 0{9^)]. Thus for a 
prolate body (A/ < 0), such as an American football, the 
body precesses slower than it spins, while for an oblate 
body the inertial precession frequency is higher than the 
spin frequency. Since the denominator in (2.5) is also 
proportional to A/, the wobble angle decreases regard- 
less of the sign of this factor. This viewpoint will be 
useful when we consider the radiation reaction problem 
for an elastic body. 



III. RADIATION REACTION FOR RIGID 
BODIES: LOCAL FORCE 

We will now re-derive the spin-down and alignment 
timescales by adding the Burke-Thorne local radiation 
reaction force to the equations of motion. 

The Burke-Thorne radiation reaction potential at a 
point X is given by m: 



B.R 



G 



rX X 



,d^Iab 
dt^ 



(3.1) 
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where lab denotes the trace-reduced quadrupole moment 
tensor: 



fab = I p{XaXb - \^^abX^) dV. 



(3.2) 



Note that this is related to the moment of inertia tensor 
according to 



lab 

with the result that 



lab - i;Sab / px dV, 
J Jv 



AI = h-Ii = -{h-ii] 



(3.3) 



(3.4) 



The radiation reaction force (on a particle of unit mass) 
is Fi^^ = —d^^^/dx"". The instantaneous {not time- 
averaged) torque on a body can easily be shown to be 



(3.5) 



Making use of Eq. ( |3.4| ) it is straightforward to calculate 
this torque for the free precessional motion. We find 

2C 

T = — A/205 sin 61(16 sin^ 6* + cos2 6l)n_Lrid , (3.6) 

acting always in the plane containing the angular mo- 
mentum and the symmetry axis X3, and perpendicular 
to Ud, i.e. along the direction of n±na shown in Fig. 1. 
We will refer to this plane as the reference plane. 



nd 




FIG. 1. For the rigid body the gravitational radiation reac- 
tion torque T lies in the reference plane. It acts perpendicular 
to the symmetry axis, i.e. along the direction of unit vector 



The evolution equations can be calculated without go- 
ing to the trouble of writing down Euler's equations. Dif- 
ferentiation of = J/Ii gives 



and so 



_ J 

^ IT' 

rsin( 
h 



(3.7) 



(3.8) 



Define J_Ln^ as the component of the angular momentum 
perpendicular to the symmetry axis. Then differentiation 
of the trivial relation 



sm( 



J 



leads to 



J 



T cos 61 



(3.9) 



(3.10) 



where Tlj is the component of the torque perpendicu- 
lar to J. Equations (B.S) and (3.10) show that the ac- 



tion of the torque breaks down neatly into two parts. 
The component along J acts to change the inertial pre- 
cession frequency while the component perpendicular 
to J acts to change 9. Substitution of (3.6) into (3.S) 
and (3.10) th en r eprodu ces the spin-down and alignment 
of equations (2^) and {2S), so the two methods of cal- 
culation agree. As t his torq ue for mulation makes clear 
(by combining Eqs. 3^ and 3.1C| ), the product 4' cos 9 
remains constant, so that if a body is set into free pre- 
cession described by {9o,<j)o): it tends to a non-precessing 
motion about X3 with (inertial frame) angular velocity 
= cos ^0 00 • 



IV. TORQUE-FREE PRECESSION OF ELASTIC 
BODIES 

We now review the theory of the free precession of an 
elastic body. This problem was first addressed in the con- 
text of the Earth's own motion. A rigorous treatment of 
the methods employed can be found in Munk and Mac- 
Donald . The terrestrial analysis was extended to neu- 
tron stars by Pines and Shaham [|. The energy loss due 
to gravitational waves was considered by Alpar and Pines 
§. 

Following the latter authors we will model a star con- 
sisting of a centrifugal bulge and a single additional de- 
formation bulge. Alpar and Pines wrote an inertia tensor 
for the elastic body of the form 

I = Iq,sS + AIninnnn - 1/35) + AldindUd - 1/35). 

(4.1) 

where S is the unit tensor [1,1,1], no is the unit vector 
along the star's angular velocity il, and is the unit 
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vector along the body's principal deformation axis (ex- 
plained below). The /o,s and A/^ pieces of I together 
represent the inertia tensor for the corresponding non- 
rotating star. The Aid term is just the non-spherical 
piece of this tensor (approximated as axisymmetric) . If 
the star were a perfect fluid, Aid would vanish, but in 
real stars (and the Earth) Aid is non-zero due to crustal 
shear stresses and magnetic fields. The term A/ji (> 
and oc fl^ for small fl), represents the increase in the 
star's moment of inertia (compared to the non-rotating 
case) due to centrifugal forces. Since the crust of a ro- 
tating NS will tend to "relax" towards its oblate shape, 
having Aid > is surely the typical case in Nature. (E.g., 
if one could slow the Earth down to zero angular veloc- 
ity without cracking its crust, it would remain somewhat 
oblate: the crust's "relaxed, zero-strain" shape is oblate, 
and after centrifugal forces are removed, the stresses that 
build up in the crust will act to push it back towards 
that relaxed shape.) But a negative Aid is also possible 
in principle. We say the deformation bulge aligned with 
Ud is 'oblate' if Aid > and 'prolate' if Aid < 0. 

What is a typical magnitude for Aid in real, spinning 
NS's? Let us assume Aid is due primarily to crustal 
shear stresses (as opposed to stresses in a hypothetical 
solid core, extremely strong B-fields, or pinned super- 
fluid vortices). Then for a relaxed crust (i.e., a crust 
whose reference ellipticity is very close to its actual ellip- 
ticity), we have Aid = bAIfi, where Alpar and Pines ||] 
estimate b ~ 10"*^ for a primordial (cold catalyzed) 
crust. The maximum value for Aid/ 1 is therefore of or- 
der 10~^. The parameter b (which arises from inter- 
nucleon Coulomb forces) scales like the average Z'^/A 
of the crustal nuclei. Since crusts of accreted matter 
(as in LMXB's) have smaller-Z nuclei [|l0|, their b fac- 
tor is correspondingly smaller, by a factor '^2 — 3. Us- 
ing AIq/I ^ 0.3(z/s/kHz)^, we would therefore estimate 
Ald/I ^ 10^^ for a NS with a relaxed, accreted crust 
and Vs 300 Hz, while for the Crab one would expect 
Aid/ 1 ~ 3 X 10^^ (again, assuming its crust is almost 
relaxed). For the freely precessing pulsar reported in 
Stairs et al. , where the body- frame precession period 
is ~ 2 X 10^ times the rotation period, Eq. (4.15) below 
(vahd for elastic bodies) yields Aid/ 1 = 5 x 10~^. For 
6 = 10~^ this corresponds to a reference oblateness of 
5 X 10"''. This is consistent with the star's crust having 
solidified when it was spinning at about 40 Hz, assuming 
that neither glitches nor plastic flow have modified its 
shape since. (When the effects of crust-core coupling are 
taken into account, giving Eq. (|7.5| ), this initial frequency 
reduces to 12 Hz. See Jones ||11[] for a review of pulsar 
free precession observations). 

Precession occurs when Ud and no are not ahgned. 
Below we describe the precessional motion when there is 
no damping. This analysis is quite general: it applies to 
any star whose inertia tensor is described by Eq. (4.1), 
independent of what causes the deformation bulge. In the 
case of several equally important sources of deformation 
along different axes, extra terms must be added to (4.1) 



and the analysis would become more complex. 



To proceed it is necessary to use equation ( 4T ) to form 
the angular momentum J of the body. However as we are 
not modelling a rigid body, we must take care to allow 
for relative motion of one part with respect to another. 
Following [Q we will write the velocity of some point in 
the body as the sum of a rotational velocity with angular 
velocity $1 and a small velocity u relative to this rotating 
frame. We will call the frame that rotates at fl the body 
frame although it is only in the rigid body limit that the 
body's shape is fixed with respect to this frame. In other 
words the velocity of some particle making up the body 
is the sum of the body frame velocity Jl x r at that point 
r plus the velocity u of the point relative to the body 
frame. Then 



Ja = lab^b + ha, 



(4.2) 



where the possibly time- varying moment of inertia is de- 
fined in the usual way: 



lab = / p{XcXcSab " XaXb) dV, 

'v 



(4.3) 



while ha is the angular momentum of the body relative 
to this frame: 



K = / peabcXbUc dV. 



(4.4) 



We will neglect the hi term when constructing a free 
precessional motion, as it can be shown that hi is small 
in a well-defined sense [|ll|. Therefore we will simply 
write 



Ja — ^ab^b 



(4.5) 



Having formulated the problem in this manner it is 
straightforward to show that the free precession of an 
elastic body is similar to that of a ri gid one. F irst write 
down the angular momentum using (4.1) and (4.5). Re- 
ferring all of our tensors to the body frame, with the 
3- axis along Ud'- 

J = (/o,5 + 2/5AIn - l/3AId)n + Am^Ud. (4.6) 



This shows that J , and Ud are coplanar. As the angu- 
lar momentum is constant this plane must rotate about 
J. As in the rigid body case, we will refer to this as the 
reference plane. See figure |[ Taking the components of 
(4.6) we obtain: 



Ji = (/o,s + 2/iAIn - l/3A/d)r!i = hVL^, (4.7) 
■h = (/o,s + 2/3A/n - l/3A/d)r!2 = /il^a, (4.8) 

.h = (lo.s + 2/3A/0 + 2/3A/d)r!3 = hQs. (4.9) 

These equations show that despite the triaxiality of I the 
angular momentum components themselves are struc- 
turally equivalent to those of a rigid symmetric top. The 
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equations of motion of the body (i.e. Euler's equations) 
involve onl y th e com ponents of J and fl. Therefore 
equations ( [4.7[ ) — ( |4.9D show that the free precession of 
the triaxial body is formally equivalent to that of a rigid 
symmetric top. We can think of the elastic body as hav- 
ing an effective moment of inertia tensor diag[Ii^Ii,I^]. 
Note that the effective ohlateness I3 — Ii is equal to A/^. 

Now introduce standard Euler angles to describe the 
body's orientation, with the polar axis along J. Let 9 
and (f) denote the polar and azimuthal coordinates of the 
deformation axis, while ^ represents a rotation about 
this axis. We refer to 9 as the wobble angle. Taking the 
ratio of components Ji and J3 using (g^ and (|^) at 
an instant when 0,2 — we obtain 



n 



bnj + iprid- 



(4.13) 



tan 7 = — tan0, 
where 7 denotes the (17, n^) angle. See figure 



(4.10) 




Oblate 



Prolate 




FIG. 2. This shows the reference plane, which contains the 
deformation axis nd, the angular velocity vector fl and the 
fixed angular momentum J. The vectors Ud and rotate 
around J at the inertial precession frequency (j>. The terms 
'oblate' and 'prolate' refer to the deformation bulge. 



We will label the angle between J and O as 9: 



(4.11) 



This angle is m uch smaller than 9, as can be seen by 
linearising ( 4.1C| ) in A/o and A/^ to give 



A/rf 

/3 



(4.12) 



Note that according to our conventions, when the defor- 
mation bulge is oblate A/^ and 9 are positive, but when 
the deformation bulge is prolate A/^ and 9 are negative. 
We can decompose the angular velocity according to 



Substituting this into equation (4^) and resolving along 
nj and gives 



Aid 

h 



(4.14) 



(4.15) 



where J denotes the magnitude of the angular momen- 
tum. Note that when A/q = the above formulae reduce 
to the familiar rigid body equations. 

Thus the motion is simple. As viewed from the in- 
ertial frame the deformation axis rotates at a rate (f> 
in a cone of half-angle 9 about the angular momentum 
vector. This angular velocity is sometimes called the 
inertial precession frequency. The centrifugal bulge ro- 
tates around the angular momentum vector also, but — 
for oblate deformations — on the opposite side of J, mak- 
ing an angle 9 = j — 9 with J. Superimposed upon this is 
a rotation about the deformation axis at a rate ip, known 
as the body frame precession frequency or sometimes sim- 
ply the precession frequency. This frequency is negative 
for an oblate distortion and positive for a prolate one. 



V. RADIATION REACTION FOR AN ELASTIC 
BODY: ENERGY AND ANGULAR MOMENTUM 
BALANCE 

Here we derive the wobble damping time tq for elas- 
tic bodies, based on energy and angular momentum bal- 
ance. Once fully underway, the derivation is just a couple 
lines. But to understand it, it is useful to carry along a 
simple, physical model for the deformed crust. (How- 
ever our derivation will actually be completely general). 
Here is the model: take some non- rotating, spherical NS, 
and stretch a rubber band around some great circle on 
the crust. We shall refer to this great circle as the NS's 
equator. Obviously the effect of the rubber band is to 
make the NS slightly prolate (but still axisymmetric) . 
To get an oblate shape, you can instead imagine sewing 
compressed springs into the surface of the crust at the 
equator. For definiteness, let the potential energy of the 
band (or springs) he V = ^eP , where I is its length. So 
e is positive for the rubber band (prolate deformation, 
A/rf < 0) and negative for the springs (oblate deforma- 
tion, A/d > 0). Now give the NS angular momentum J 
about some axis that is not quite perpendicular to the 
equator. We now have our deformed, wobbling NS. We 
consider the equation of state of the star and the value e 
to be fixed once and for all, and consider how the energy 
of the system (star -I- band) varies as a function of its 
total angular momentum J and the wobble angle 9 (the 
angle between J and the perpendicular to the equator); 
i.e., we consider i?( J, 9). We will be concerned with small 
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wobble angle, so let us expand E{J,9) as a Taylor series 
in J and 9: 



E{Ja 



1 



1 

24' 



1 



(5.1) 



Here Eq is defined to be the energy of the (star + band) at 
zero J, and _B, C, and F are some expansion coefficients 
that in principle depend on the physical properties of the 
(star + band). Fortunately we will soon see that there 
are simple relations between B, C, and F and previously- 
defined physical parameters, such as A/^. Our ultimate 
goal is to obtain the tw o pa rtial derivatives on the right- 
hand side of equation ( |2.5| ), where E now denotes the 
total energy. 

First, to see that no lower order terms (such as 
J, 6>J , 0^ , or O.P terms) can appear in the expansion 
( ^.iD , note that the J — configuration corresponds to 
the minimum of the potential energy of the (star -I- band) 
system. Displacements of the (star -I- band) are first or- 
der in .P, so changes in the potential energy of (star -|- 
band) are 0{J^). Thus terms in E{J,9) that are oc 
are kinetic energy pieces. These terms with a in them 
are clearly just ^{lQ^)'^''JaJb, where /g*' is defined to be 
the inertia tensor of the (star -I- band) at J = 0. (Cor- 
rections to the star's 1°"^ first enter the energy at ©(J"*).) 
We write /q** as 



TO:b J cab I A T f ^ah ^ cab 



(5.2) 



where /o,5 represents the 'spherical part' of /q . Then 



1 



cab 



Aid 

lo.s 



(5.3) 



where a term of 0{AI^) has been neglected. The kinetic 
energy part of E is [up to terms of 0{AI^) and 0( J^) ] 



J2 



210,3 



Ah 

Io,s 



(5.4) 



where we have used the small wobble angle result Ja^"^ 



J(l — h9^). From Eq. (5.4) we i mmed iately read off the 



values of B and Fe in expansion (5.1) 



2 Aid. 
Fe = AId/{Io,sf 



and obtain the partial derivative 

Aid 



dE 



(5.5) 



(5.6) 



To compute the partial derivative in the numerator of 
equation ( ^.5| ) it is sufficient to consider the 6 ^ Q limit 
i| so that 



dE 

in 



}=0 



BJ+-CJ^, 
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(5.7) 



where f2 denotes the spin frequency in the axisymmetric 
limit. It is related to the inertial precession frequency by 



n^^ii-Aid/io.s)- 

The final physics inputs we need are 
2G A/2 



E 



5c5 /o„ 



E = (t)J . 



(5.8) 



(5.9) 



(5.10) 



Eqs. (5.9) and ( pi.lC| ) follows from the quadrupole for- 
malism in the same way as for the rigid body. 

The necessa ry p ieces have been gathered; substituting 
into equation (2^) gives 



ej^ Aid 

2G A/, 



(5.11) 
(5.12) 



This is simply the same spin-down rate as for a rigid 
body, with the replacement (A///i) — > ed- This is much 
longer than the timescale claimed by Bertotti and Anile 
1^ by a factor Ala/AId, which is typically ~ 10^ or 
higher. 

Finally, the spin-down rate can be obtained in the 
same way as for a rigid b ody, i.e. b y diff erentiating (f) — 
J/Ii and using equations (5.9) and (5.10). Strictly there 
will also be a term in Ii , but this correction will be down 
by a factor of order {n/flmax)'^- We then obtain the same 
spin-down as for a rigid body, again with the replacement 
AI ^ Aid-. 



Io,s 



(5.13) 



VI. 



RADIATION REACTION FOR AN ELASTIC 
BODY: LOCAL FORCE 



We now give a second derivation of the wobble damp- 
ing rate for an elastic star, by directly adding the gw 
radiation reaction force to the Newtonian equations of 
motion. Besides being a satisfying consistency check on 
the calculation in §4, by doing this second derivation cor- 
rectly we can show where Bertotti and Anile ||^ went 
astray. 

As was the case for the rigid body, the Burke-Thorne 
potential will exert a torque on the spinning star. How- 
ever, this is not the only effect of the radiation reaction 
force: It will distort the shape of the NS and thus its mo- 
ment of inertia. The equation describing the precession 
is then of the from 



-[ilN + SlBT)m=T, 



(6.1) 
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where In denotes the Newtonian part of the moment 
of inertia tensor, SIbt the perturbation in this tensor 
due to the Burke-Thorne force, and T the Burke- Thorne 
torque. It was the SIbt terms that were not included 
by Bertotti and Anile. Fortunately, these can also be 
calculated explicitly, as we show below. 



the unit vector in the reference plane which lies perpen- 
dicular to J and points towards n^. We then find 



RR 



5c5 



AldO - MnO]{vanjh + njaVb) ■ 

(6.4) 



A. Effect of on the NS's Shape 

It is perhaps surprising that one can explicitly deter- 
mine the effect of on the NS's moment of inertia, 
since the answer would seem to depend on the NS's mass 
and the details of its equation of state; i.e., one might 
worry that extra parameters must be specified even to 
make the problem well-defined. However the point is 
that (from symmetry arguments) the perturbation A/jj 
depends only on a single physical parameter, and this 
parameter already appears in our Newtonian equations 
of motion. That parameter is A/o/f2^, the amount of 
oblateness caused "per unit centrifugal force" . 

The point is that both the centrifugal and radiation 
reaction forces have the very special property that they 
grow linearly with distance from the center of the star. 
This fact, coupled with symmetry arguments, is enough 
to determine A/y in terms of A/n/f^^; no new physical 
parameters have to be introduced. 

Let $^ be some external potential of the form (^'^ = 
A°'^XaXb, where Aa& is some trace-free tensor. Allow this 
potential to act on the non-rotating (and so spherically 
symmetric) NS; it will induce a perturbation A/"'' in 
the NS's inertia tensor. Since the background is spher- 
ically symmetric, the only possibility (to first order in 
the perturbation) is that A/"'' = CA"'', where C is some 
constant (i.e, independent of A°^). 

We can determine C as follows. Decompose the cen- 
trifugal potential into a spherically symmetric and a 
trace-free piece: 



1 



1 



VL^x^ + K'^XaXb, (6.2) 



where A^" 



For small Vl the per- 



turbed inertia tensor is A/"'' — AIii{nf^n'^ — ■|(5°''), so 
the constant C is just 2A/o/il^. 

The radiation reaction potential for the freely precess- 
ing elastic body can be found by substituting the radia- 
tion reaction free motion into equation (3.1) to give: 



RR 



X X 



^Id-^{ndandb) + AIn — {nnannb) 

(6.3) 



The first term is the potential caused by the motion of the 
deformation bulge, the second by the centrifugal bulge. 
The differentiations of the unit vectors are straightfor- 
ward. In the case where <^ 1 we can approximate 
nd ~ nj + dn±j and n^i « nj — 6n±j, where n±j is 



Here v denotes a unit vector nj x n±j. Using the pre- 
scription described above, these radiation reaction po- 
tentials can be converted immediately into perturbations 
of the moment of inertia tensor: 

2C 

SIbt = -^^^[A/rfA/o^ - iAIn)^e]{vnj + njv) . 

(6.5) 



B. Adding <I'^^ to Equations of Motion 

It now remains to compute the torque T using equa- 
tion (3^). We obtain four terms, corresponding to the 
expansion of the product of I with its fifth time deriva- 
tive. Again linearising with respect to 9 we obtain 



2C 



e + AIdAIne~Alle]n^j. 

(6.6) 



Define = A/o/Iq.s and e<j = Ald/Io.s- Then the 
terms on the rhs of ( |6.6| ) stand in the ratio ea/^n ■ '■ 
1 : eo . We are now in a position to write down the equa- 
tion for d{I N^)/dt. Using Eq. ( |6.5| ) and the Newtonian 
motion to compute d\{SlBT) ^]/d t, and neglecting terms 
of order 9^, we find that Eq. (|6.l|) reduces to 



l(/.n) 



-^^'0[AhAIne 



{AInfd]n^j (6.7) 



2G : 
5? 



^(j,^[Alje - AldAlJ + AldAInO - All^§]n±j. 



We see that the last two terms on the rhs are cancelled 
by terms on the Ihs. This leaves 



dt 



2G_ 

5c5' 



[Alle - AIdAlJ]n^j. (6.8) 



The problem has reduced to a rigid-body Newtonian one, 
with the two torque terms indicated on the right-hand 
side. The terms stand in the ratio 1 : en. In fact, the 
dominant term is the same as that obtained in the rigid 
body case with the change AI Aid- 



^ Note our definition of en differs by a factor 2/3 from 
who set en = |A/n//o,s- 
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We therefore find that the ahgnment rate as calculated 
using the local Burke-Thorne formalism agrees with the 
flux-at-infinity method. The previous force-based cal- 
culation of Bcrtotti and Anile |^ failed to include the 
deformation 6Ibt, so that the cancellations in equation 



( p. 71) described above did not occur. 

Finally, it is easy to show that even when the approx- 
imations <IC 1, Ed <C 1 are not employed, the effective 
torques due to the 6Ibt terms are still perpendicular 
to J, so the spin-down (j) using this local formalism is 
necessarily the same as in the flux-at-infinity method. 



and 



dE 
'do 



'61 A/. 



(7.3) 



(See Jones for a detailed derivation.) 

These lead to an alignment timescale that is /crust 
shorter than that of equation ( ^.13 ). This result is con- 
firmed using the local torque formulation, where 



T 



(7.4) 



VII. ALLOWANCE FOR A LIQUID CORE 

We have successfully described the effects of gravita- 
tional radiation reaction on an elastic precessing body. 
We will now briefiy describe how to extend this result to 
the realistic case where the star consists of an elastic shell 
(the crust) containing a liquid core. The Earth itself is 
just such a body, and the form of its free precession was 
considered long ago. We will base our treatment on that 
of Lamb , who considered a rigid shell containing an 
incompressible liquid of uniform density. To make the 
problem tractable the motion of the fluid was taken to 
be one of uniform vorticity. We will assume the elliptic- 
ity of the shell, and also the ellipticity of the cavity in 
which the fluid resides, are small. Then the small an- 
gle free precession of the combined system can be found 
by means of a normal mode analysis of the equations of 
motion |12|. 

The key points are as follows: The fluid's angular ve- 
locity vector does not significantly participate in the free 
precession. Instead it remains pointing along the sys- 
tem's total angular momentum vector. The shell pre- 
cesses about this axis in a cone of constant half-angle. 
The fluid exerts a force on the shell such that the shell's 
body frame precession frequency is increased in magni- 
tude, so that: 



■0 = -0- 



A/ 



(7.1) 



where A/ denotes the difference between the 1 and 3 
principal moments of inertia of the whole body, not just 
the shell. 

We now wish to calculate the alignment rate of such a 
body due to gravitational radiation reaction. The aver- 
aged energy and angular momentum fluxes, as well as the 
instantaneous torque, depend only upon the orientation 
of the mass quadrupole of the body, and so are ex actl y 
t he same as i f the body were rigid, i.e. equations (2.1), 
(2.2) and (3.6) apply. Equations giving the kinetic energy 
and angular momentum of the body are given in Lamb 
. These can be used to obtain the partial derivatives 
that appear in equation (ETS). Explicitly, we find 



dE 

a7 



vt = 4) + 4> 



(7.2) 



In the realistic case where both crustal elasticity and 
core fluidity are taken into account we can combine the 
above arguments as described by Smith and Dahlen , 
i.e. we can take the rigid result and put / /crust and 
A/ Ala to give 



~(f> 
2G A/2 



1 A/d 



5c5/, 



crust 



(7.5) 
(7.6) 



VIII. CONCLUSIONS 

We have shown that the gw damping time for wobble 
in realistic NS's has the same form as for rigid bodies. 



but with the replacement AI^/Ii A/^//c 
given an alignment timescale of: 



This 



1.8 X 10^ yr 



/„ 



10^4 g cm^ 



10'^**gcm^ 
A/rf 



/kHz 



(8.1) 



For the Crab, taking ^ 3 x 10~^, this gives Tg ^ 
5 X lO-'^^yrs-much longer than the age of the universe. 
For an accreting NS with td ^ 10^^ and Vg ~ 300 Hz, we 
estimate Tg ~ 2 x 10® yrs. 

Our basic conclusion, then, is that gw backreaction is 
sufficiently weak that other sources of dissipation proba- 
bly dominate. Unfortunately, even for the Earth the dis- 
sipation mechanisms are not well understood Q. Early 
estimates of Chau and Henriksen [Q, which considered 
dissipation within the neutron star crust, suggested that 
wobble would be damped in around 10^ free precession 
periods, i.e. over a time interval of 10^/(edJ^s). A more 
recent study of Alpar and Sauls ||l^ argued that the dom- 
inant dissipation mechanism will be due to imperfect cou- 
pling between the crust and the superfluid core. They 
estimate that the free precession will be damped in (at 
most) 10^ free precession periods. In contrast, according 
to equation (Bj), the gw damping time is in excess of 

10® ^Mis^ free precession periods. On the basis of these 

estimates, it seems likely that internal damping will dom- 
inate over gravitational radiation reaction in all neutron 



9 



stars of interest. Note however, that while internal dissi- 
pation damps wobble for oblate deformations, we expect 
that internal dissipation causes the wobble angle to in- 
crease in the prolate {Aid < 0) case. 

A study of the gravitational wave detectability of re- 
alistic neutron stars undergoing free precession, includ- 
ing a discussion of other astrophysical mechanisms which 
might affect the evolution of the motion, will be presented 
elsewhere (Jones, Schutz and Andersson, in preparation). 
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